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Thermal quantum and classical correlations in two qubit XX model in a nonuniform 

external magnetic field. 

Ali Saif M. Hassan, 1 -!] Behzad Lari, 2 ^Q and Pramod S. Joag 2 -| 

1 Department of Physics, University of Amran, Amran, Yemen 

k Department of Physics, University of Pune, Pune, India-4 11007. 

(Dated: December 9, 2010) 

We investigate how thermal quantum discord (QD) and classical correlations (CC) of a two 
qubit one-dimensional XX Heisenberg chain in thermal equilibrium depend on temperature of the 
bath as well as on nonuniform external magnetic fields applied to two qubits and varied separately. 
We show that the behaviour of QD differs in many unexpected ways from thermal entanglement 
(EOF). For the nonuniform case, (Bi = — B2) we find that QD and CC are equal for all values of 
(Bi = —B2) and for different temperatures. We show that, in this case, the thermal states of the 
system belong to a class of mixed states and satisfy certain conditions under which QD and CC 
are equal. The specification of this class and the corresponding conditions are completely general 
and apply to any quantum system in a state in this class and satisfying these conditions. We 
further find that the relative contributions of QD and CC can be controlled easily by changing the 
relative magnitudes of B\ and i?2. Finally, we connect our results with the monogamy relations 
between the EOF, classical correlations and the quantum discord of two qubits and the environment. 
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It is now quite well known that composite quantum 
systems can be in a class of states, called entangled states, 
in which the correlations between the constituents of the 
system cannot be achieved in a classical world (l|, l2| ■ Al- 
though all pure entangled states possess such nonlocal 
quantum correlations, there are mixed entangled states 
which do not, in the sense of violating Bell inequalities 
[3| . The entanglement in quantum states and the result- 
ing nonlocal quantum correlations form an area of in- 
tense research, due to their huge technological promise, 
especially in the areas of quantum communication and 
cryptography [J]. However, quantum correlations break- 
ing Bell inequalities need not account for all quantum 
correlations in a composite quantum system in a given 
state. In order to account for the quantum correlation 
in a given state, we must find some means to divide the 
total correlation into a classical part and a purely quan- 
tum part. This is particularly important for mixed states, 
since their quantum correlations are many a time hidden 
by their classical correlations (CC). An answer to this 
requirement is given by quantum discord (QD), a mea- 
sure of the quantumness of correlations introduced in Ref. 
5] . Quantum discord is built on the fact that two classi- 
cally equivalent ways of defining the mutual information 
turn out to be inequivalent in the quantum domain. In 
addition to its conceptual role, some recent results |6j, 
suggest that quantum discord and not entanglement may 
be responsible for the efficiency of a mixed state based 
quantum computer. Present authors believe that QD will 
turn out to be a very useful tool to analyze mixed state 
quantum correlations and their consequences, as mixed 
state entanglement is very difficult and eluding to deal 
with [7|. To realize this hope we need a viable relation 
between mixed state entanglement and quantum discord 



8]. The pointers towards such a relation may be ob- 
tained by studying these properties for various quantum 
systems. 

Motivated by these considerations, we present here the 
results of our investigation of the amount of QD and CC 
in a two qubit Heisenberg XX chain at finite temper- 
ature subjected to nonuniform external magnetic fields 
B\ and B2 acting separately on each qubit. We study 
two distinct cases namely Bi = —B 2 (nonuniform field) 
and Bi = B 2 (uniform field). In each case, we obtain the 
dependance of QD, CC and entanglement of formation 
(EOF) [l7| in the system on the external magnetic field 
and temperature. Such a model is realized, for example, 
by a pair of qubits (spin 1/2) within a solid at finite tem- 
perature experiencing a spatially varying magnetic field. 
Such Heisenberg models can describe fairly well the mag- 
netic properties of real solids [9j and are well adapted to 
the study of the interplay of disorder and entanglement 
as well as of entanglement and quantum phase transi- 
tions JlO, 11 1 . The variation of entanglement [12| and 



QD [13|, [14| in a two qubit Heisnberg XX chain with 



external magnetic field is already reported. 

In order to quantify entanglement in a thermally mixed 
two qubit state, we use entanglement of formation de- 
rived using concurrence given by [15] 



C = max{Ai — A 2 — A 3 — A 4 , 0} 



(1) 



where Xi(i — 1, 2, 3, 4) are the square roots of the eigen- 
values of the operator pp in descending order 



p=(oX®ol)p*(oX®ol), 



(2) 



with Ai > A2 > A3 > A4, and p is the density matrix 
of the pair qubits; a\ and a\ are the normal Pauli op- 
erators. The entanglement of formation is related to the 



concurrence by 



EN = h 



i + VT~c^ 



where h(x) — — a;log 2 x— (1 — x) log 2 (l — x). Henceforth, 
in this paper, we denote the entanglement of formation 
(EOF) by EN. The concurrence C — corresponds 
to an unentangled state and C = 1 corresponds to a 
maximally entangled state. 

The thermalized Heisenberg system. The model Hamil- 
tonian we study is given by 



H = J(SfS% + SfSl) + B^l + B 2 S Z 2 



(3) 



where S a = a a /2, (a = x, y, z) are the spin 1/2 opera- 
tors, cr a are the Pauli operators and J is the strength of 
Heisenberg interaction. B\ and B 2 are external magnetic 
fields. As stated in the introduction, by changing B\ and 
B 2 separately, we want to study the effects of magnetic 
field on the thermal QD, CC and EN in a very general 
way. The eigenvalues and eigenvectors of H are 



H\QQ) = -(B 1 +B 2 )\00) ) 



£r|ll) = (Bi + B 2 )|ll) 



H\^)=±D\^), 



(4) 



where D 2 = (B x - B 2 f + J 2 and |^±) = ^-[|01) 



( 1- , 1 10)]. We denote the eigenvalues correspond- 
ing to |00), |X1), |-0 =1= ) by Eoo,En,E± respectively. In the 
standard basis, {|00), |01), |10), |11)}, the density matrix 
p(T) is given by 



P{T) = 



Z 



where 



U\ 



mi 

mi v 

v w 2 

u 2 



,(B 1 +B 2 )/kT 



(5) 



u 2 = e 



-(Bi+B 3 )/kT 



and Z = Tr[exp(-rff-)] is the partition function. In the 
following we select \J\ as the energy unit and set k = 1. 

Quantum Discord, p, Il3| In classical information theory 
(CIT) the total correlation between two systems (two 
sets of random variables) A and B described by a joint 
distribution probability p(A, B) is given by the mutual 
information (MI), 



I(A, B) = H(A) + H(B) - H(A, B) 



(7) 



with the Shannon entropy H(-) = —J2jPj^°9^Pj- Here 
Pj represents the probability of an event j associated to 
systems A, B, or to the joint system AB. Using Bayes's 
rule we may write MI as 



I(A,B)=H(A)-H(A\B), 



(8) 



where H(A\B) is the classical conditional entropy. In 
CIT these two expressions are equivalent but in the quan- 
tum domain this is no longer true [5|, [l6| . The first quan- 
tum extension of MI, denoted by I(p), is obtained directly 
replacing the Shannon entropy in Eq.(7) with the von 
Neumann entropy, S(p) = —Tr(plog 2 p), with p, a density 
matrix, replacing probability distributions. To obtain a 
quantum version of Eq. (8) it is necessary to generalize 
the classical conditional entropy. This is done recogniz- 
ing H(A\B) as a measure of our ignorance about system 
A after we make a set of measurements on B. When B 
is a quantum system the choice of measurements deter- 
mines the amount of information we can extract from it. 
We restrict ourselves to von Neumann measurements on 
B described by a complete set of orthogonal projectors, 
IIj, corresponding to outcomes j. 

After a measurement, the quantum state p changes to 
Pj = [(I ®Ylj)p(I ®Hj))/pj, with / the identity operator 
for system A and pj = Tr[(I ® Wj)p(I ® n,)]. Thus, one 
defines the quantum analog of the conditional entropy as 
S^pl-jTlj}) = ^2jPjS(pj) and, consequently, the second 
quantum extension of the classical MI as 5] 

J(p\{nj}) = S{p A ) - S{p\{n 3 }). The value of 
J{p\{Rj}) depends on the choice of {H,}. 

Henderson and Vedral [5( have shown that the max- 
imum of i7(p|{H,-}) with respect to {H,} can be inter- 
preted as a measure of classical correlations. Therefore, 
the difference between the total correlations I(p) and 
the classical correlations Q(p) — supsji^J(p\{Hj}) is 
defined as 



, . D . (B 1 -B 2 ) , , D , 
Wl =cosh(-)+ - smh(-), 



D (gi - B 2 ) . D 

W2 =cosh(-) — smh (-), 



Jsinh(T§.) 



D(p) = I(p) - Q(p), 



(9) 



D 



giving, finally, a measure of quantum correlations [5| 
called quantum discord (QD). For pure states QD re- 
duces to entropy of entanglement 17f , highlighting that 
in this case all correlations come from entanglement. 
However, it is possible to find separable (not-entangled) 
mixed states with nonzero QD [5|, [l6| , meaning that en- 
W tanglement does not cause all nonclassical correlations 



contained in a composite quantum system. Also, QD 
can be operationally seen as the difference of work that 
can be extracted from a heat bath using a bipartite sys- 
tem acting either globally or only locally jig . 
Results and discussion. Case I: B\ = —B2, and (J > 0). 

In this case, \ip~) is the ground state with eigen- 
value E- = —-\/AB'l + J 2 . Other eigenvalues are 0.4 
{0, 0, y/4B'( + J 2 } for eigenvectors {|00), |11), \ip+)} re- 
spectively. 

In this case the variation of QD, CC and EN with B\ 
at different temperatures (T=0.2, 0.9, 1.5) is depicted 
in Figs. 1, 2, 3, respectively. We observe that QD and 
CC in the thermal state coincide for all values of B\ for 
different temperatures (T=0.2, 0.9, 1.5). 

In order to understand this observation, we take a close 
look at the thermal state. At temperature T, the thermal 
state is given by 



1 



P=-[|00)<00| + |11)(11| 



^IBI + JZ/T 



|^~) vTI + e-V 4B ?+^ a /r^+)^+|] , 



This p has the Bloch representation [21 [ 



P=l\I< 



y^ Cjo-j 



(10) 



(11) 



where e^ (i = 1,2,3) are the one-qubit Pauli opera- 
tors. In the appendix to this paper we prove that, the 
class of mixed states as in Eq.(ll) have equal classical 
and quantum correlations (like bipartite pure states [19() 
provided 



C.j > Cfc 



and 



('k 



(12) 



where i / j / I; £ {1,2,3}. Here Ci,Cj,Ck are the 
diagonal elements of the correlation matrix defined by 
Cij = Tr(p<Ji ® o-j). It is straightforward to check that 
the thermal state p in Eq.(10) which has form of Eq.(ll) 
satisfies conditions (12). This explains the observations 
in Figs. 1, 2, 3, that the two qubit thermal state p for 
B\ = —B 2 in Eq.(10) gives rise to equal QD and CC for 
all values of B\ and temperature. In order to see why 
the common curve for QD and CC peaks at B\ = we 
can maximize the expression for QD = CC with respect 
to B\ and check that the maximum occurs at B\ = 0. 

From Figs. 1, 2, 3, we also see that EN as a function 
of B\ has a peak at B\ = for T = 0.2, has a dip for 
T = 0.9 and goes to zero over an interval symmetric 
about B 1 = for T = 1.5 [12|. From Eq. (10) we see 
that the concurrence of the thermal state is governed by 
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Figure 1 : QD and CC (dashed line) and EN (solid line) as 
a function of external magnetic field B\ = — B2 at T = 0.2 




Figure 2 : QD and CC (dashed line) and EN (solid line) as 
a function of external magnetic field Bi = — B2 at T = 0.9 




Figure 3 : QD and CC (dashed line) and EN (solid line) as 
a function of external magnetic field Bi = — B2 at T = f .5 



the admixture of the \ip + ) and l^ - ) states. We expect 
concurrence to fall as the state \ip + ) classically mixes 
more and more with the ground state \tp~)- For fixed J = 
1 and a fixed temperature T, this happens for B\ = 0. 
That is why we got a dip in the EN curve at B\ — 0. The 
size of this dip increases with temperature. In fact the 
dip touches the B\ axis when, at B\ = the concurrence 
is zero. To see this we note that the concurrence for p in 



Eq.(lO) is given by J2(j 



C = — maxlH - y/uiu 2 ,0} 



(13) 



where v, u\, u 2 are given in Eq. (6). Therefore, for B\ = 
and J = 1, C > provided sinh ± > 1 or T < 1.1346. 
For T = 1.5 and J = 1, using the requirement sinh ^ = 
13, we can find the range of B\ around B\ = in which 
C = 0. This is -1.1456 < Si < 1.1456. Figs. 1,2,3 
confirm the corresponding behaviour of EN. 

Fig. 4 shows the variation of EN, QD, CC with temper- 
ature at fixed values of B\ = —B 2 . As expected we have 
QD = CC for all temperatures. Both EN and QD (CC) 
curves have plateau at low temperatures corresponding 
to their ground state values, as at these temperatures, the 
ground state is not thermally connected to other exited 
states. Other interesting observation is the vanishing of 
concurrence at a finite critical temperature T c which in- 
creases with Bx value, while QD and CC asymptotically 
go to zero with temperature. The increase in T c with B\ 
[l2| can be understood from the thermal state Eq.(10) 
which says that higher temperatures are required to get 
a given admixture of \ip~) and \ip + ) for higher B\ values. 




Figure 4 : QD and CC (dashed line) and EN (solid line) as 

a function of the absolute temperature T for (a) 

Si = -B 2 = 1 (b) B 1 = -B 2 = 2 



We now deal with the case B 2 = —aB\, a ^ 1 and pos- 
itive. B\ — satisfies both, B 2 = —aB\ and B 2 = —B\, 
so that QD = CC at B\ = 0, for all temperatures T. For 



a fixed temperature T, it turns out that QD > CC for 
Si ^ if a > 1 and CC > QD for B x ^ if < a < 1. 
This is depicted in Fig. 5, for a = 2 and a = 1/2 for 
T = 1.5. The dominance of QD over CC (or vice versa) 
varies continuously with a. This observation gives us the 
key to control the contributions of QD and CC to a two 
qubit thermal state in Heisenberg model via the continu- 
ous variation of the applied magnetic field. The behavior 
of concurrence in this case can be analyzed in a way sim- 
ilar to the case B 2 = —Bi,(a = 1). From Fig. 5, we 
see that for the same temperature, the range over which 
concurrence vanishes depends on a, this range decreases 
monotonically with a. Also, the peak position of con- 
currence (or, EN) on the B\ axis shifts monotonically 
towards B\ — as a increases. Thus the main entan- 
glement features can be controlled by varying external 
magnetic fields. 
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Figure 5 : QD (dashed line) , CC (dash-dotted line) and EN 

(solid line) as a function of external magnetic field Si at 

T = 1.5 (a) B 2 = -2Bi (b) B 2 = -Bi/2. 



Case II : B x = B 2 . 

For uniform external magnetic field B\ = B 2 , Figs. 6, 
7, 8, show the variation of QD, CC and EN with B\ for 
temperatures T = 0.2, 0.9, 1.5, respectively. We see that 
all three quantities are symmetric about B\ = where 
they have their maxima. Further, QD > CC, except at 
B\ = 0, where QD = CC. For higher temperatures, the 
qualitative behavior of QD and CC remains the same, 
while EN curve drops down below those of QD and CC. 
This can be qualitatively understood by looking at the 



thermal state given by 
1 



P 



^[e 2Sl / T |00)(00| 



,-2Si/Ti 



11)(11| 



+e J / 7 >-)W>-| + e - J / 7 >+)(V; + |]. 



(14) 



For small temperatures, the entanglement of the ther- 
mal state is largely dictated by that of l^ - ) and becomes 
dominant. At higher temperatures, admixture due to 
other states reduces the entanglement, so that QD and 
CC dominate. Such a complementary behavior of en- 
tanglement and discord can serve as a pointer towards a 
possible connection between them. 




Figure 8 : QD (dashed line) , CC (dash-dotted line) and EN 

(solid line) as a function of external magnetic field Bi where 

Bi = B 2 at T = 1.5 
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Figure 6 : QD (dashed line) , CC (dash-dotted line) and EN 

(solid line) as a function of external magnetic field B\ where 

Si = B 2 at T = 0.2 
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Figure 7: QD (dashed line), CC (dash-dotted line) and EN 

(solid line) as a function of external magnetic field B\ where 

Si = B 2 at T = 0.9 



Fig. 9 shows the variation of QD, CC and EN with 
temperature for B\ = B^ = 1 and B\ = B^ = 2. We 
see that for high temperatures, QD hugely dominates 
EN, showing the robustness of QD with temperature. 
As temperature becomes large QD and CC converge 
towards each other. For larger values of B\ this happens 
at higher temperatures. As the temperature increases, 



all the coefficients in the thermal mixture Eq.(14) tend 
to be equal and the thermal state approaches random 
mixture. Thus it seems that QD and CC approach 
each other as an arbitrary thermal state approaches 
a random mixture. Obviously, for random mixture 
p = \ (I I), QD = CC = 0. A quantitative analysis of 
the relative behaviors of QD and CC with temperature 
will be very interesting, but possibly have to wait for 
further developments in the theory. 
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Figure 9 : QD (dashed line) , CC (dash-dotted line) and EN 

(solid line) as a function of the absolute temperature T for 

(a) Bi = B 2 = 1 (b) Si = B 2 = 2 
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Figure 10 : ENab (solid line), QD^ (dashed line), EN AE 

(dot line) and QD<^ (dash-dotted line), as a function of 

external magnetic field Si = — B2 at T — 0.9 



It will be intresting to connect our results with the 
monogomy relations between the EN and the classical 
correlation 22] of two subsystems (qubits) and the envi- 
ronment 



EN A1 



cc a~e - Sa > 



EN AE + CC^ = S A , 



and the relation between EN and QD [22 



EN AB + EN AE = QD 



:w 



QD 



AE 1 



(15) 



(16) 



showing us that EN and QD always exist is pairs. Here 
A,B label the qubits and E stands for the environment. 
We assume that environment (heat bath) comprises the 
universe minus the qubits A and B so that the state 
Pabe is a pure state. Since the variation of all the 
quantities pertaining to the system AB with B\ and T 
are obtained form the XX model, we can use Eq. (15,16) 
to find the corresponding dependence of EN A e and 
QD^, on B x and T. Figs. 10 and 11, (for B x = -B 2 ) 
show the variation of EN AE and QD<t-=, with B\ and T. 



The monogamic relations also help us establish a nec- 
essary and sufficient condition for QD*^ = CC^ when 
the environment is present. This is : QD^ — CC+rs if 
and only if \l AB = EN AE + EN AB - QD^. To prove 
the necessity we note that when QD^ = CC*r-g, (that 
is, QD-^-g = \l AE ), Eq.(16) can be written as 

-Iab - EN AE + EN AB - QD^. (17) 

Now suppose Eq.(17) is true. Then using Eq.(16) we have 

Q D A-B = \^B 



which implies QD<r-=, = CC*r-^. Figs. 12 and 13, show the 
variation of both sides of Eq.(17) with B\ and T which 
establishes Eq.(17) for the XX model. 
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Figure 11 : ENab (solid line), QD-^ (dashed line), ENae 

(dot line) and QDvg (dash-dotted line), as a function of the 

absolute temperature T for B\ = —B2 = 1. 
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Figure 12 : \l AB (solid line), ENae + ENab - QD^ 

(dashed line), as a function of external magnetic field 

B 2 = -B\ at T = 1.5 




Figure 13 : \Iab (solid line), ENae + ENab - Q-D^g 

(dashed line), as a function of the temperature T for 

B 2 = -B x = -1. 



Summary. In this paper we have studied the variation 
of QD, CC and EN in two qubit XX Heisenberg chain 



as functions of independently varied magnetic fields 
B\ and B 2 on each qubit and also with temperature. 
We deal with two cases Si = — B 2 (nonuniform field) 
and Bx = B 2 (uniform field). Our first observation 
is the complementary behavior of entanglement and 
QD/CC ' . For the nonuniform magnetic field, we get 
the interesting observation that QD and CC are equal 
for all B\ = —B2 values as well as all temperatures. 
Surprisingly, this observation is explained quite simply, 
using the symmetric form of the thermal state. A very 
interesting observation is that the relative contributions 
of QD and CC can be tunably controlled by varying 
the applied magnetic field. Another interesting finding 
is that the equality oi QD and CC of the subsystem 
(qubits) imposes a constraint on the distribution of 
QD and EN over the subsystem and its environment. 
Further investigation of general Heisenberg models like 
XXZ along these lines may turn out to be interesting 
and fruitful. 

Appendix. We prove the following statement. 

If the quantum state has the Bloch representation [21| 



H [/< 



3 



(Al) 



and Ci = Cj > Ck and Ck = —cf where i 7^ j 7^ k G 
{1,2,3}, then this state contains the same amount of 
quantum and classical correlation (QD — CC). 

Proof: In Ref. 24] S. Luo evaluated analytically the 
quantum discord for a large family of two-qubit states, 
which have the maximally mixed marginal and their 
Bloch representation is 



-[I®I- 



3 

E 

i=l 



CjO-j 



For this class of quantum states the quantum mutual 
information is 

Z(p) = i[(l-ci-c 2 -C3)/o52(l-ci-C2-c 3 ) + (l-ci + 
c 2 + c 3 )log 2 (l - ci +C2 + C3) + (l + ci-c 2 + c 3 )log 2 (l + 
c\ - c 2 + c 3 ) + (1 + c\ + c 2 - c 3 )log 2 (l + ci+c 2 - c 3 )]. 

We substitute the conditions above in the quantum 
mutual information. Puting c = c\ = c 2 > c 3 and c 3 = 
— c 2 , we get, 

Ap) = |[(1 - 2c + c 2 )Zo 52 (l - 2c + c 2 ) + (1 - c 2 )Zo 92 (l - 
c 2 ) + (1 - c 2 )/o. 92 (l - c 2 ) + (1 + 2c + c 2 )log 2 (l + 2c+ c 2 )] 

After some algebraic simplification, we get 
l(p) = (1 - c)log 2 (l -c) + (1 + c)log 2 (l + c) 
which equals ICC as in Ref. .24]. It is also easy to check 
that the above argument goes through when c = C\ = 
C3 > c 2 and c 2 = — c 2 and when c = c 3 = c 2 > c\ and 
cj = -c 2 , to get l(p) = ICC. Thus, 

QD(p)=l(p)-CC(p) = CC(p). 
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